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We consider the particle creation process associated with a quantum field x in a time-dependent, 
homogeneous and isotropic, classical background. It is shown that the field square x^, the energy 
density and the pressure of the created particles have large fluctuations comparable to their vacuum 
expectation values. Possible eflFects of these fluctuations on the reheating process after inflation 
are discussed. After determining the correlation length of the fluctuations in two different models, 
corresponding to the decay in the parametric resonance regime and in the perturbation theory, it 
is found that these fluctuations should be taken into account in the flnal thermalization process, 
in the back-reaction effects and when the formation of primordial black holes is considered. In 
both models, by comparing quantum and thermal fluctuations with each other it is observed that 
very quick thermalization after the complete inflaton decay is not always possible even when the 
interaction rates are large. On the other hand, when the back-reaction effects are included during 
the preheating stage, the coherence of the inflaton oscillations is shown to be lost because of the 
fluctuations in . Finally, we note that a large fluctuation in the energy density may cause a black 
hole to form and we determine the fraction of total energy density that goes into such primordial 
black holes in the model of preheating we consider. 



I. INTRODUCTION 

The theory of quantum fields in curved spacetime is 
by now a well established subject (see e.g. [f-3]). The 
issues like the non-uniqueness of the vacuum and the par- 
ticle creation process, the renormalization of the stress- 
energy-momentum tensor are all thoroughly understood, 
and the current research is mainly focused on the incor- 
poration of interactions and alternative rigorous defini- 
tions of the theory (see e.g. [4-6]). There are also some 
unsolved problems, the primary example being the black 
hole information paradox, but it seems that these can 
be explained in the full quantum theory of gravity and 
not in the semi-classical approximation which treats the 
background classically. 

Recently, we have pointed out in [7, 8] an alternative 
description of the cosmological particle creation process, 
which tries to overcome two potential issues which may 
arise in the standard treatments based on the Fourier de- 
composition. The first one is a (possible) conflict between 
causality, which requires the particle creation process to 
take place independently in each Hubble volume, and the 
introduction of the globally defined Fourier modes as the 
main physical observable. The second issue is related to a 
subtlety in the interpretation of the vacuum expectation 
value of the number operator of a momentum mode, since 
one can show that the number of created particles with a 
fixed momentum largely fluctuates about its mean value. 
To resolve these two issues, in [7, 8] we have analyzed the 
particle creation process using wave packets localized in 
a Hubble volume. As a result, the issue of causality is 
naturally settled by treating the creation of the modes 
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in different Hubble volumes independently and the mean 
values are interpreted as statistical averages over distinct 
Hubble volumes. By applying this formalism to reheating 
process after inflation, in [7, 8] it is shown that there ex- 
ists small density perturbations on Hubble length scales 
at the end of reheating. Although the same result can 
also be obtained in the standard formulation by intro- 
ducing a suitable window function to probe the horizon 
scale, the construction presented in [7, 8] is physically 
more direct and transparent. 

Since field theory respects locality by construction, the 
first issue mentioned above should not be a problem as 
long as the modes are physically interpreted properly. 
Indeed, the fact that field variables must commute at 
space-like separations, i.e. [x{t' ,x'),x{'^-i^)] — when 
g^iyAx^ /S.x^ > 0, guarantees causality. The Fourier 
mode expansion of the field variable x(t, x) is just a con- 
venient way of describing physics. Any other complete 
set can be used for expansion and physics should not de- 
pend on the basis chosen. However, the second matter 
related to the interpretation of expectation values can be 
tricky, especially if there are large fluctuations about the 
mean values. 

In this paper, we consider the particle creation effects 
in a time-dependent, homogeneous and isotropic, clas- 
sical background and try to determine the magnitude of 
the fluctuations of the energy-momentum tensor T^i, and 
the field strength x^ about their mean values. Namely, by 
defining the fluctuation operators 5T^^ = T^i/— < T^^, > 
and 5^2 — < >, whose expectation values van- 
ish < (5T^^ >— < 5^2 >= 0, we calculate the variances 
< {ST^^Y >7^ and < {5^2^ >^ 0. As we wiU see 
in section II, for most of the components of the energy- 
momentum tensor including the energy density and the 
pressure, and for x^ the deviations have the same order 
of magnitude as the corresponding mean values, which 
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shows the existence of large fluctuations. We emphasize 
that the spatial scale of these fluctuations is given by the 
correlation length of the quantum field X- 

In the second part of this paper, in section III, we dis- 
cuss possible implications of these results for the reheat- 
ing process in single scalar field inflationary models. We 
consider the decay of the inflaton field (j) into bosonic x 
particles in two different prototype models correspond- 
ing to perturbativc and parametric resonance regimes, 
and in both models we determine the correlation length 
of the X field excited by the inflaton oscillations. We find 
that the existence of these fluctuations may affect various 
events during reheating. It is pointed out in section III 
that thermalization process cannot be completed unless 
quantum fluctuations become smaller than the thermal 
fluctuations in equilibrium and this may delay the final 
moment of reaching thermal equilibrium. Secondly, in 
the model with preheating, when the x field starts affect- 
ing the frequency of the background inflaton oscillations, 
the fluctuations in x^ ^-re found to destroy the coherence 
of the oscillations. We show that because of this loss of 
coherence it is no longer possible to neglect the spatial 
derivatives of the inflaton when the back-reaction effects 
become important. Finally in section III, the possibility 
of black hole formation because of the fluctuations in the 
energy density of created particles is considered. For a 
large fluctuation, which of course happens rarely, the cor- 
responding Schwarzschild radius may become larger than 
the correlation length which would produce a black hole. 
We show that the observational constraint on the frac- 
tional energy density that can go into primordial black 
holes puts some new restrictions, especially in models of 
preheating. We finally conclude by reviewing our findings 
and discussing future directions in section IV. 



II. THE FLUCTUATIONS 

We consider a real scalar field x propagating in a cos- 
mological Robertson- Walker metric 

ds^ = -dt^ + a^{dx'^ + dy^ + dz^), (1) 

which has the following action 

5 = -^/v^[(Vxf + MV]- (2) 

We assume that in addition to the scale factor a, the 
mass parameter M may also depend on time: M = M{t). 
Therefore, particle creation can be induced both by the 
evolution of the metric and by the externally varying time 
dependent mass parameter M. The energy- momentum 
tensor can be obtained from (2) as 

T^^u = V^xV.x - \9^.. [V/3xV^x + M\^] . (3) 

More explicitly, one can read the energy density p, the 
pressure P, the "momentum" Ui and (trace-free) stress 



components Tj,- as 
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Ui = Tu = xidiX) = \ {mx) + (SiX)x] , 

Tij = 9jX5jX- [/'(5feX)(9ix)] , 



(4) 



where dot denotes time derivative and all indices refer 
to the obvious coordinate basis of (1). Note that the 
spatial components of the energy-momentum tensor are 
decomposed as T!^- = Tjj -|- Pgij, where cj^^Tjj =0. In 
quantum theory, there is an ordering ambiguity in Ui 
and in (4) we perform a symmetric ordering. 

For quantization it is convenient to define a new field 
Xby 

X = a^l\. (5) 
Then, the action (2) up to surface terms becomes 



X' - g'^{diX){djX) - [M' - ^-H' 
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where H = a /a is the Hubble parameter. The momen- 
tum variable 11 conjugate to X is given by 11 = JsT and one 
can easily apply the standard canonical quantization pro- 
cedure. Introducing the mode function X^ and time in- 
dependent ladder operator corresponding to a fixed time 
to, ak, as 



X 



1 



(27r)3/2 



d-'k 



(6) 



the canonical commutation relation [X{t,x),'n{t,x')] = 
i5{x — x') can be satisfied by imposing 



[ak,4']='^(k-k'), 



and 



(7) 



(8) 



where boldface letters k and x refer to the spatial 3- 
vectors ki and a;', k.x = kiX^ and k"^ = 6^^kikj. On the 
other hand, the field equations imply 



where 



Xk + i^kXk = 0, 



(9) 



(10) 



Because of the homogeniety and isotropy of the back- 
ground, the mode function Xk depends only on the mag- 
nitude k of k. 
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The ground state of the system at time to can be de- 
fined as 



ak|0 >= 0, 



(11) 



and in terms of the mode functions this corresponds to 
the initial conditions 



^fe(*o) 
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I, Xk{to) 



(12) 



Using (12), the Hamiltonian at time to can be found as 



H{to) = / d-'k 



1 



(13) 



which justifies the identification of |0 > defined in (11) 
as the ground state. Let us note that the above formu- 
lation is identical to the standard formulation of particle 

creation in terms of Bogoliubov transformations (sec [9]). 
The correspondence can be achieved by defining the time 
dependent and /3k coefficients as 
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i f ujkdt' ^ if ujhdt' 
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By imposing the following equations 

Wk 2i f* ujkdt 

Oik = 7, — e ■''o 

■ Ldk -2i r* ujkdt 

Pk = 7, — e -"o ak, 

and the initial conditions ak{to) = 1, (3k{to) = 0, the 
equivalence of both formulations can easily be estab- 
lished. For future use let us also calculate Xk as 



Xu = 



—i r* utkdt' if* Ukdt' 

—lak e ■'*o -|- ipk e -"o 



(15) 



Note that because of the equations obeyed by and 
in taking time derivative of X^ from (14) one can treat 
ttk and pk as if they are constant parameters. 

The particle creation process is usually described by 
giving the vacuum expectation value of the number op- 
erator for a mode of momentum k. Alternatively, one 
can also calculate the vacuum expectation value of the 
energy-momentum tensor < T^i, >=< 0|r^j/|0 >. One 
benefit of working with the energy-momentum tensor is 
that the potential problem mentioned in the introduc- 
tion, namely the creation of a global Fourier mode in an 
expanding universe, is naturally evaded. Using (6) in (4), 
one can straightforwardly calculate 



<p> = T + V + G, 
<P> = T-F-^G, 



(16) 



<Ui> 



0, 



0, 



where being the real functions of time T, V and G are 
given by 



T = 



V = 



G = 



-r^-^ ( d^k\Xk-'^HXk?, 
2(27ra)3 J ' 2 ' 
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d?kM^\Xk^, 



(17) 
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Here, T and V play the roles of the kinetic and potential 
energies, respectively, and G is like the energy stored in 
the gradient of the field. It is an easy exercise to show 
that when M is constant the expectation value of the 
energy-momentum tensor is conserved 



M{t) = Mo 



< Tt"" >= 0, 



(18) 



which actually reduces to a single nontrivial equation 

< p> + 3H{< p > + < P >) = 0. Thus, in that case 

< T^y > can be fed into the right hand side of the Ein- 
stein's equations to include the back-reaction effects. 

The functions T, V and G defined in (17) diverge in 
general. Therefore, they should be regularized to be in- 
terpreted physically. Adiabatic (or WKB) regularization 
is a suitable and physically appealing way of getting fi- 
nite expectation values for < >, < p > and < P > 
[10, 11]. One can see that the function V is related to 

< > and thus it can be made finite by adiabatic reg- 
ularization. Similarly, the functions T and G can also be 
expressed in terms of < p >, < P > and < x^ >> thus 
they can also be regularized. When we mention about 
the magnitudes of the variables T, V and G below, we 
implicitly refer to their regularized finite values. 

The fact that < Ui >= and < Tij >= does not 
imply Ui = and Tij = identically. Similarly, p and P 
do not exactly equal their mean values. All these physical 
quantities fluctuate about their expectation values and 
we would like to find out the size of these fluctuations. 
For that we deflne the fluctuation operator 

ST^. = T^u- < T^. >, (19) 

and try to compute the variance < ((5T^iy)^ >^ 0. 

Let us start with the energy density dp = p— < p >. 
Using (6) in the definition of p given in (4) and reading 

< p > from (16), a relatively long but straightforward 
calculation gives 



<6p'^> = 2T2-K 21/2+ ^2 

o 



(20) 



-I- 



M2 



(27ro)6 

M2 



d^kcos{cl>k)\Xk- -HXk\\Xk\ 



d^ksm{(f>k)\Xk--HXk\\Xk\ 



(27ra)6 
the ph; 

and Xk- The only "nontrivial" step in this calculation 



where (j>k is the phase difference between {Xk — ^HXk) 
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is to deal with a double integral which has the following 
form 



(21) 



This can be converted by suitable angular integrations 
into 



(22) 



which can then be expressed in terms of the function G. 
As a check on this computation one can see from (17) 
and (20) that 



KSp^ > > 2T2 + 2y2 

M2 



d'k\Xk- -HXk\\Xk\ 



(23) 



(27ra)6 

= n/o^ ^fi / d^ki(fk2[ak,ak2-bk,bk2f>0, 
2[2Tra)° J 

where Uk = \Xk — ^HXk\ and bk = M\Xk\. This shows 
that < Sp^ >> 0, as it should be. 

Let us now turn to the fluctuations of the pressure 
and define SP = P— < P >. Again a straightforward 
calculation, which uses (4), (6) and (16), gives 



< <5p2 > = 2T^ + 2V^ + — 

27 



(24) 



+ 



(27ra)6 
(27ra)6 



/ 



d'kcos{cPk)\Xk-^HXk\\Xk\ 



d^ksm{<f>k)\Xk- -HXk\\Xk\ 



Here also one has < 6P^ >> 0, as expected. 

Both in (20) and in (24), there are integrals involv- 
ing the phase <j)k, which cannot be expressed in terms 
of the functions T, V and G. Interestingly, however, 
these terms appear with opposite signs in < Sp^ > and 
< (5P2 >, and from (20) and (24) one sees that 

20 

<5p'^> + < 5P^ >= 4T2 + 4y2 + _g2_ (25) 

Comparing with the average < p > given in (16), eq. 
(25) shows that < 5p'^ > + < SP^ > has the same order 
of magnitude as the mean energy density square 



<dp'^> + < 6P'^ >= c(< p >f 



(26) 



where c is a number close to unity, which varies depend- 
ing on the hierarchy between the functions T, V and G. 
For instance, c ~ 32/27 if T ~ G > or c ~ 0.97 if 
T ~ F - G. 

Before proceeding with the calculations of other fluc- 
tuations, let us make a few comments. It is easy to see 

from (23) that the integrals involving cos((/)fe) and sm{(j}}S) 
in (20) and (24) are bounded by 2T'^ + 2V'^. Therefore, 



one would expect them to be rcgiilarizcd in a suitable 
scheme. On the other hand, the sum of the fluctua- 
tions < 5p^ > + < 5P^ > can already be expressed 
in terms of finite functions T, V and G after adiabatic 
regularization. In general, one would also expect the os- 
cillating cos((/)fc) and sin(^fe) integrals in (20) and (24) 
to be smaller than the previous terms. Therefore, both 
< 6p'^ > and < SP"^ > should have the same order of 
magnitude as (< p >)2. In any case, (25) shows that 
the fluctuations < Sp^ > and < SP^ > cannot be simul- 
taneously small compared to the average energy density. 
This can be viewed as an uncertainty relation between 
the energy density and pressure fluctuations. 

Because of these large fluctuations, it may not always 
be possible to use a simple equation of state P = wp to 
characterize the energy-momentum tensor in quantum 
particle creation. One can calculate the correlation of 
the energy and pressure fluctuations as 



< {5p, 5P} >= 4T2 - 4y2 



\G' 



where the symmetric ordering {5p, 5P} — 6pSP + SP5p is 
chosen inside the brackets. Thus, depending on the case 
the correlation can be negligible, which would forbid the 
use of an efltective equation of state parameter w. One 
can also determine the expectation value of the commu- 
tator, which should give information about simultaneous 
measurability of fluctuations. We find 



< [SP,Sp] >-- 



4M2 



J 



'^kcos{cPk)\Xk--HXk\\Xk\ 



d^ksin{^k)\Xk- -HXk\\Xk\ 



therefore the integrals involving the phase 4'k is related 
to the expectation value of the commutator. Note that 
the expectation value is equal to zero for M = 0. 

Let us now continue with the other components of the 
energy- momentum tensor, namely Ui and Tij. As noted 
in (16), their expectation values vanish. To determine 
fluctuations, one can calculate 



< UiUj >= -TGgij. 



and 



|G2 



ifG2 if i = j, 



(27) 



(28) 



where in the last equation the summation convention is 
not used, i.e. i and j are treated as free indices. The cor- 
relation of field variables carrying different set of indices 
can be seen to vanish because of the homogeneity and 
isotropy of the background. However, when G has the 
same order of magnitude as T and V, the fluctuations 
of the stress components r^j cannot simply be ignored. 
Remember that G is the energy stored in the gradient 
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of the field, therefore it is not surprising to discover that 
the fluctuations of the stress components depend on G. 
Similarly, for large G and T, the "momentum" compo- 
nents Ui has large fluctuations, of the order of average 
energy density. Again it is natural to see both T and 
G in the fluctuations of momentum, since T measures 
kinetic energy and G measures spatial variations. 

When back-reaction or symmetry breaking effects are 
considered, becomes an important parameter and 
thus it is crucial to determine the fluctuations about the 
mean < >■ Defining as before the variation operator 
S^2 = x^— < x^ > a-iid using (5) and (6), one can easily 
calculate 



(29) 



Therefore, viewing x^ ^ random variable one sees 
that the corresponding standard deviation is equal to \/2 
times the average, which again indicates the existence of 
large fiuctuations about the mean. This result can indeed 
be anticipated without doing any computation since the 
free field x(^jx) can be viewed as a collection Gaussian 
random variables defined at each point in space and (29) 
is true for any Gaussian distribution. 

Let us remind that the expectation values of physical 
quantities determined above should be interpreted as sta- 
tistical averages over different points in space at a given 
time. However, the field variables at nearby points are 
not independent, i.e. they are correlated with each other. 
The size of a region containing such correlated field vari- 
ables is given by the correlation length ^c, which can be 
determined from the two point function 



<x(i,x)x(t,0) >= 
1 



1 



27r2a3 



(27ra)3 
dk k^lXkf 



d^k\Xk\^e'^-'' 
sin(fcr) 



kr 



(30) 



where = SijX^x^ . The comoving correlation length 
is the minimum value of r that makes (30) to vanish. 
Physically, it gives the spatial comoving size of a region 
in which the x field is homogenous. Similarly, it is also 
possible to define the correlation lengths of the energy 
density, the pressure etc., however these are expected to 
have the same order of magnitude as ^c- Therefore, one 
can imagine the space to be divided into comoving re- 
gions of typical size ^c, and in each region the physical 
quantities like the x field, the energy density p and the 
pressure P are uniformly distributed. The fluctuations, 
on the other hand, give variations from one region to 
another (see figure 1). 

As r — >■ 0, (30) should give the (regularized) mean 
value < x^ >. Therefore, in calculating the correlation 
length from (30), a suitable regularization should be per- 
formed before the momentum integral is taken (note that 
(30) is finite for r 7^ without any need of regulariza- 
tion). To employ WKB regularization [11] one uses (14) 
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FIG. 1: The space divided into regions of typical (comoving) 
size ^c. At a given instant, the x field can be assumed to 
vary appreciably from one region to another, while it is nearly 
uniform in each region. 



in (30) and sets ja^p = -f 1, which yields 



<x(i,x)x(i,0) >= 



dk — X 



(31) 



1/3.1^+ Re {a.Ple-'' 



sin(fcr) 
kr 



where an additive factor of 1/2 is ignored in the first 
square brackets for regularization. Since initially one im- 
poses Pkito) — 0, the momentum integral in (31) equals 
zero at i = to and it is expected to converge also at later 
times. We will use this regularization in determining the 
correlation length in the next section. 

Let us note that in our calculations we ignore the 
back-reaction effects and assume that the evolution of 
the background does not change as the particles are cre- 
ated. This is an approximation and, for instance, as the 
energy of the created particles increases the energy of 
the background should decrease. Naively, one may then 
tend to view these fiuctuations as " isocurvature" pertur- 
bations, however, this depends on how the energy density 
of the background and the created particles redshift. On 
the other hand, keeping in mind that the back-reaction 
effects are ignored may help to solve the following con- 
cern. As noted in (18), when M is constant the average of 
the energy- momentum tensor is conserved and one may 
worry that this can be spoiled if fluctuations are added 
on top of the mean values. However, all that is needed 
is to satisfy the energy-momentum conservation for the 
whole system, and this should be guaranteed when the 
back- reaction effects are properly taken into account. 

It is interesting to compare the direct evaluation of the 
dispersion used in this paper with the standard formula- 
tion which employs a smearing of the dispersion (or power 
spectrum) in momentum space with a suitable window 
function to probe a given scale. In general these two 
computations should agree, at least on large scales, how- 
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ever there arc some important differences. While using a 
window function can be a convenient way of regularizing 
some expressions, it is not possible to identify the corre- 
lation length, since the power spectrum does not contain 
any information about the correlation of variables at dif- 
ferent points. Moreover, some results may be sensitive 
on the choice of a window function whereas the adia- 
batic regularization offers both technically and physically 
unique way of getting finite and consistent values from 
singular expressions. To emphasize the difference, let us 
point out that recently in [30, 31] adiabatic regularization 
is applied to cosmological perturbation theory and some 
standard results, such as the form of the power spectrum, 
are shown to be modified significantly depending on the 
regularization. 



III. IMPLICATIONS FOR REHEATING 

In this section, we analyze the particle creation process 
in the reheating period after infiation and consider the 
single scalar field inflationary models. In general, the 

evolution of the metric and the inflaton </) is governed by 
the Friedmann and the scalar field equations 



3M| 



■ dV 
3i^0+-=O, 



(32) 



where V{(f>) is the scalar potential. In reheating, the 
scalar oscillates about the minimum of the potential and 
we focus on models in which the potential at that stage 
can be taken as 



V = -m 
2 



2 J,2 



(33) 



where m is the infiaton mass. Because of the expansion 
of the universe, the amplitude of the oscillations slowly 
decreases in time so one can assume a solution of the 

form 



(j) = <^{t) sm{mt). 
When 6 <C m$, the field equations (32) become 



-H^ = 0, 



and these can be solved as 

2/3 



$ = 



Swrnt 



(34) 



(35) 



(36) 



different types of interactions suitable for preheating and 
perturbative decay, respectively. We first determine the 
correlation length of the excited x field at the end of 
decay and evaluate the functions T, V and G introduced 
in the previous section in (17). Later we study possible 
implications of our findings on the three important pro- 
cesses, which are the final thermalization of the decay 
products, the back-reaction effects and the formation of 
primordial black holes. 



Preheating 

Let us first start with preheating. The decay of the 

inflaton fleld in the parametric resonance regime has been 
analyzed in detail in [12-15]. Specifically, it is shown in 
[13, 15] that the following interaction 



Hi 



2/2 2 

;g <P X , 



(37) 



can give a decay due to broad parametric resonance if 
the parameter q =^ g^^^/rn^ ^ 1. This condition can 
naturally be satisfied in chaotic inflationary scenario in 
which $0 ~ Mp, where $o is the initial value of the 
inflaton amplitude and Mp is the Planck mass. In this 
model the frequency (10) becomes 



^1 



fc2 



+ 5^$^ sin^(TOi), 



(38) 



and the particle creation mainly occurs due to the oscil- 
lating term. As shown in [13, 15], for momenta in certain 
instability bands, the solution of (9) exponentially grows, 
which corresponds to jafcj ~ \/3k\ ^ 1- Actually in the 
expanding universe the real process is much more compli- 
cated where a given momentum, which is initially in the 
first band, jumps over many different instability bands in 
time. Fortunately, it can still be treated analytically by 
approximating the process in terms of successive scatter- 
ing on parabolic potentials [15]. The final result for the 
Bogoliubov coefficients can effectively be described by an 
index fj,/. such that 



(39) 



Moreover, the first and the most important resonance 
band is initially peaked around the physical momentum 
fc* = x/gm^O) with a width of the order of fc* (here we 
only consider the first stage of preheating where back- 
reaction and re-scattering effects are ignored). The effec- 
tive index Hk can be approximated as 



(40) 



Thus, the evolution is equivalent to the dust dominated 
universe and in that case the combination 9if^/4-|-3if/2, 
which appears in (10), is equal to zero. 

We assume that reheating occurs due to the coupling 
of the infiaton to a bosonic field, x, and consider two 



where /ij,(fc,) ~ 2/i/(/c*)^, /i is the average index depend- 
ing on the coupling constant g and it typically varies 
in-between 0.1 and 0.2 [15]. 

To calculate the correlation length we use (39) in (31). 
As discussed in [15], the integral of the second term in 
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the square bracket in (31) (the one containing the real 
part) gives a time dependent oscillating correction which 
is less than 1 compared to |/3fe|^ integral (see (89) in [15]). 
Therefore, (31) can be written as 



<x(i,x)x(i,0) >c 



27r2a3 



dk 



sm(kr) 
kr 



(41) 

where C is the correction factor mentioned above. Using 
(39) and (40), the momentum integral in the last equa- 
tion can now be performed using the steepest-decent ap- 
proximation, which gives 



<x(i,x)x(t,0) >c 



sin(fcH.r) 



k^T 



<X > ■ 



(42) 



Therefore, the comoving correlation length can be deter- 
mined as 



1 



1 



(43) 



In this computation, the effects of the modes in other 
instability bands and the time evolution of the first band 
caused both by the redshift of momenta and the decrease 
of the inflaton amplitude which would lessen the width 
of the band, are ignored (note that fc* given below (39) 
is the physical momentum). However, as shown in [15], 
the modes that have been amplified from the very begin- 
ning become exponentially larger than the others. Thus, 
in practice one can treat fc* as a comoving momentum 
scale such that the first instability band does not change 
in time and (43) should give a good estimate for the co- 
moving correlation length of fluctuations. 

It is interesting to compare the size of the correla- 
tion length to the Hubble radius at the end of the pre- 
heating. The broad parametric resonance ends when 
q = g^f^"^ /m? ~ 1 and this can be used to estimate the 
value of the amplitude at that time as $ ~ m/g. The 
Hubble parameter can be determined from (35) as 



H 



m 



(44) 



By comparing the initial and the final values of q, the 
amplitude $ can be seen to decrease ^/qo times, where 
go is the initial value of q at the beginning of preheating. 

Then, from (36), one sees that the universe expands q^^^ 
times and thus the physical correlation length at the end 
of broad parametric resonance is given by 



phys 



1/3, 



1/3 
% 



1/12 _i 

= % m . 



(45) 



As a result, (43) and (44) gives the ratio of the physical 
correlation length to the Hubble radius as 



^phys 



R 



H 



5/12 ^ 



(46) 



One can have $o ^ but (jo ^ 1 in the broad para- 
metric resonance regime, thus the ratio is in general less 
than unity. 

In this model, the variables T, V and G can be deter- 
mined in terms of < x^ > as follows. Firstly, from (17), 
one sees that 



1 



< X >= 2^ * sin^(mi) < x > 



(47) 



Therefore, V oscillates between zero and the maximum 
value Vmax — g^^^ < > /2. To determine G, one can 
again use (39) and (40) in (17), and apply the steepest- 
decent approximation. The steps are identical to the 
derivation of (42) and one finds 



G 



k^ 



< X > 



(48) 



To find T, one can first ignore the expansion of the uni- 
verse to a very good approximation. Using (15) in (17) 
one then obtains 



T 



1 



Re akPte J*o 



2i r* oji^dt' 

tn 



dk 



where we regularize this expression by ignoring an addi- 
tive factor of 1/2 in the square brackets. Applying again 
the steepest-decent approximation, we find 



1-C 



< X >, 



(49) 



where C is the factor defined in (41). Since from (38) 
(jJk, is an oscillating function of time, T is also oscillating. 
However, it is easy to sec that V + G ^ T . 

Finally, it is also possible to determine the phase 0^, 
which first appeared in (20). Ignoring the expansion of 
the universe, (j}^. equals the difference between the argu- 
ments of Xk and Xk- Using \ak \ — |/3fe|, one can see that 
Arg(Xfc) ~ ipk/2 and Arg(Xfc) ~ ip,^/2 + tt, where ipk is 
the phase difi^erence between ak and /3fe. Therefore 



<f)k — 



(50) 



From (20) and (24), the value of (pk can be seen to fa- 
vor the energy fiuctuations compared to pressure fiuctu- 
ations. 



Perturbative decay 

Let us now consider the reheating process in perturba- 
tion theory. For that we assume the following trilinear 
coupling 



1 



O-0X 



(51) 



which may arise after spontaneous symmetry breaking. 
In this model, the frequency (10) becomes 



k^ 

ljI = + (7$ sin{mt). 



(52) 
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As discussed in [17], the perturbation theory is appUca- 
ble when a'^/m? <^ 1, which can in general be satisfied if 
the ampHtude $ is smaU. It is known that the preheating 
picture completely changes when both interactions (37) 
and (51) present in the Lagrangian, see [16]. Therefore, 
the perturbative dec;ay due to (51) should be considered 
on its own as a diffcsrcuit model, i.e. it is not to be pre- 
ceded by the preheating considered above. Let us remind 
that the evolution of the background fields is still given 
by (36). 

In perturbation theory i.e. for ]^fe| ^ 1, Bogoliubov 

coefficients can iteratively be solved and to first order 
they can be determined as 



1, 



The time integral in (53) can be evaluated using the sta- 
tionary phase method [18] (note that in almost all models 
of inflation <C m, which is important for the applica- 
bility of the stationary phase approximation) . There are 
two oscillatory terms in (53), one is the explicit pure 
phase exponential and the other is coming from Uk/iOk. 
It can be shown that for a given k, the main contri- 
bution to the integral comes from an interval near 
fixed by 0Jk{t*) = m/2, which, in the perturbative regime 
cr$/m^ <C 1, implies 



a. 



m 
2"' 



(54) 



This can be interpreted as the decay of the inflaton at 
time to two x particles with comoving momentum k 
(see e.g. [17]). Thus, only the modes in the following 
interval 



aom , aim 

-^r- < k < —— 



(55) 



significantly produced, where Uq and oi are the scale fac- 
tors at the beginning and at the end of the decay, re- 
spectively. We quote from [7] the result of the relatively 
straightforward calculation: 



/2 V. 



(56) 



2rn5/2 \2kJ ' 

where $o is the initial value of the amplitude. 

To determine the correlation length, we use (53) in (31) 
and find 



<x(i,x)x(t,0) >- 



1 



27r2a3 







dk — X 



*dt'^cos(2 J^^ Mt")dt" 



(57) 

sin(A;r) 
kr 



The time integral in the square brackets in (57) can also 
be performed using stationary phase method. However, 
the phase integral now produces an extra cosine term, 

cos I / [2^^Jk{t')+m]dt'\ , (58) 



where is the time corresponding (54). This term can 
be seen to oscillate very rapidly in the decay range (55) 
and thus the second term in the square brackets in (57) 
can be neglected since its contribution will be very small 
after performing the momentum integral. Using (56) in 
(57) then gives 



raim/2 g^iQ^ikr) 

<x(t,x)x(t,0) / dk-^, 

Jaom/2 ' 



(59) 



where B is a fc-independent constant. The indefinite in- 
tegral can be explicitly evaluated in terms of the Fresnel 
cosine integral C{x) as 




2 sin(fcr) 
Vk 



(60) 



One can now see that to a very good accuracy the comov- 
ing correlation length is independent of the upper limit of 
the integral aim/2, instead it is fixed by the lower limit 
as 



1 



aom 



(61) 



Thus, the physical correlation length at the end of the 
decay is given by 



^phys ^ 



aom' 



(62) 



which is larger than 1/m by the expansion factor of the 
universe during reheating. 

By noting from (56) that the spectrum is given by 
|/3fc| ~ k~^/^, it is easy to understand (62) physically 
as follows. The dependence of |/3fe| on the comoving mo- 
mentum k indicates that, as far as the correlation length 
is concerned, the particle creation effects decrease with 
increasing momentum. Since the process can be thought 
as the decay of the inflaton particle with mass to into two 
X particles with physical momentum to/2, the particles 
created in the beginning of the process have the small- 
est comoving momentum aom/2 and thus are the most 
important ones. The corresponding wavelength, which 
redshifts in time, gives the correlation length (62). All 
the particles created in due time have larger momenta 
and smaller wavelengths, and they produce sub-leading 
corrections to the correlation length. 

To compare the correlation length to the Hubble ra- 
dius, let us determine the Hubble parameter at the end 
of the decay. As discussed, e.g. in [17], the decay process 
described by the interaction (51) is equivalent to decay 
with a constant decay rate ~ ct^/to, which gives a re- 
heating temperature Tr ~ ^^T^Mp ~ a'^M^jm. The 
Hubble parameter corresponding to this temperature is 
given by 



a 
m 



(63) 
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Therefore, 



R 



■H 



(64) 



In these models one usually imposes a m for pertur- 
bation theory to be applicable. On the other hand, the 
expansion factor during reheating depends on the initial 
value of the inflaton amplitude, but it is not expected to 
be a very large number. Thus, one again finds that the 
correlation length is smaller than the Hubble radius. 

In perturbation theory, the functions (17) can also be 
approximately determined in terms of < >• By defi- 
nition the potential energy V is given by 



out that a very quick or instant thermalization of the de- 
cay products may also be prohibited by the existence of 
quantum fluctuations in the energy density. 

It is known that in thermal equilibrium at temperature 
T the energy density is given by 



(68) 



where 7 is a constant depending on the number of bosonic 
and fermionic species in equilibrium. Then, Uy = pV, 
where Uy is the average energy inside a physical vol- 
ume V. In the canonical ensemble, the dispersion of the 
energy about its mean value is given by 



V = -fT$sin(mt) < > • 



(65) 



The calculations of the variables T and G proceed as which, using (68), implies (sec e.g. [27]) 
follows. Compared to the momentum integral in < x^ >, 



the integrand in T and G contains two more powers of 
k (note that in the decay range ui^ — k/a). Evaluating 
the integrals, the dominant contribution comes from the 
Ultra- Violet (UV) end of the limit, which is equal to m, 
and this gives 



G 



(66) 



Since in perturbation theory (T$/m^ ^ 1, one finds that 



T ~ G » 14, 



where VmaT, is the maximum value of 



the potential energy V . 

From (14) and (15) the phase can also be deter- 
mined easily when the expansion of the universe is ne- 
glected. Using \(Xk\ ^ |/3fc|, one finds that 0^ ~ 37r/2. 
As oppose to the preheating, the integrals involving 
in (20) and (24) now tend to decrease 5p and increase 
bP. Finally, by using (56) in (31) (and taking the r ^ 
limit) one can obtain 



(67) 



Not surprisingly, one sees from (67) that in perturbation 
theory < x^ > is smaller by several orders of magnitude 
compared to Mp, and it is also much smaller than rv?. 



Thermalization 

After determining the correlation length of quantum 
fluctuations, let us now discuss possible effects that they 
might produce during reheating. We first consider the 
thermalization process of the decay products. This is 
a difficult process to study and unfortunately there is 
not much work done in the literature (see e.g. [19-26]). 
In general, the distribution of particles produced during 
the decay of the inflaton can be seen to be far from ther- 
mal equilibrium [19]. Thus, depending on the interaction 
rates the full thermalization can be delayed resulting a 
low reheat temperature. Here, we would like to point 



(5C/y _ 2 1 
W " (7^)1/2 ^W^- 



(69) 



(70) 



In (70), the dependence of the relative deviation on V is 
a characteristic feature of thermal fluctuations. 

Thermal equilibriiim can only be justified if thermal 
fluctuations at a given scale is larger than the quantum 
fluctuations at the same scale. Otherwise, one should 
wait for some time for the energy to be transferred in 
space for uniformization. (see the discussion following 
(76) below). In a volume V ^ there are on the average 
F/(^P''^*)'^ number of uncorrelated regions as far as quan- 
tum fluctuations are concerned. Since in each region the 
relative deviation of the energy density is of the order of 
unity, the relative order of quantum fluctuations in the 
volume V can be found as 



Uy 



3/2 



(71) 



Comparing with (70), one sees that instant thermal equi- 
librium after the decay can only be justified if 



T < 



(72) 



i.e. when the thermal correlation length is larger than 

^phys ^ Eq. (72) places an upper bound on the reheating 
temperatiire in terms of the correlation length of quan- 
tum fluctuations. Note that since thermal fluctuations 
can occur locally, the volume V about which the fluctu- 
ations are compared should be taken inside the horizon 
and thus it is important to have ^p^v^ to be less than the 
Hubble radius. 

Since the Hubble parameter should not change 
abruptly, the final equilibrium temperature Tr can be 
determined from the Friedmann equation 



M2 



Mr 



(73) 
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which gives Tr ~ yjHMp. The full thcrmalization can be 
achieved when ~ H , where is the total interaction 
rate of the decay products and should not be confused 
with r^, i.e. the decay rate of the inflaton. 

In the preheating model reviewed above, the Hubble 
parameter at the end of the broad parametric resonance 
is determined in (44). If thermal equilibrium sets in a 
very short time following the end of the broad parametric 
resonance, the reheating temperature must be fixed as 




Then, (45) and (72) give 



^ 2/3 1/6 

9>t' % ■ 



(74) 



(75) 



Since q^y^ 1, the coupling constant becomes large g > 1, 
which shows that the quantum fields </> and x a^re strongly 
interacting. This is a very difficult regime to study and 
the results about the preheating period reviewed above 
can no longer be trusted. For (/ < 1, (72) cannot be sat- 
isfied, which shows that even when the total interaction 
rate is large to yield quick thermal equilibrium in prin- 
ciple, this cannot be achieved since would be thermal 
fluctuations are (much) smaller than quantum fluctua- 
tions. 

On the other hand, in the perturbative decay de- 
scribed by the interaction (51), the reheating temper- 
ature is given by Tr ~ ^a'^Mpjm (here one assumes 
~ F^ ~ o-^/m). Then, (62) and (72) imply 



(76) 



It is not very difficult to satisfy this condition since one 
usually assumes ct ^ m <C Mp. 

If (72) is not satisfied, one should wait for the energy to 
be redistributed for thermal equilibrium. This process is 
different than the redistribution of mean energy or mean 
number density by scattering or particle decays, as ana- 
lyzed, for instance, in [22]. Namely, fluctuations should 
propagate in space-time to achieve uniformization and 
this cannot happen instantaneously because of causality. 

To illustrate what would happen when (72) is violated, 
i.e. when quantum fluctuations are larger than thermal 
fluctuations, assume that F^ is large enough to convert x 
particles into radiation in local thermal equilibrium in a 
very short time. Consider now two different regions that 
have the size of the correlation length and assume that in 
one region all the energy stored in the inflaton oscillations 
is converted into x particles and in the other region only 
half of the energy is converted into x particles. This is 
a typical situation since the deviation in the energy den- 
sity of the created x particles is comparable to its mean 
value. Let us now try to imagine what would happen as 
the universe expands twice. In the first region the en- 
ergy density, being composed of radiation in equilibrium. 



will be redshifted by 16 times, while in the second region 
half of the energy density, i.e. the radiation, will be red- 
shifted by 16 times and the other half corresponding to 
inflaton oscillations will be redshifted 8 times. Therefore, 
there appears a genuine density contrast 6p and one can 
easily see that because of the different redshifts of the 
constituent energy densities, Sp/p grows with the expan- 
sion like 6p/p ^ a. 

On the other hand, the free propagating of radiation 
is expected to work for the uniformization of energy, but 
this process is slowed down by causality. Namely, by 
looking at the evolution of radial null geodesies, a co- 
moving spherical region of volume V can be seen to ex- 
pand (in the comoving grid) as a^/^. Therefore, Sp/p is 
expected to decrease like Sp/p ^ 1/VV ^ a"^/"* which 
is smaller than the increase Sp/p ^ a noted above. As 
a result, one sees that the fluctuations in the x field in- 
duce changes in the total energy density during the decay 
process. Later on, when everything is converted into ra- 
diation, the fluctuations tend to be smoothed out, but de- 
pending both on the magnitude and the scale, one should 
wait for some time for full thermal equilibrium to set in. 

This naive scenario should be sharpened by consider- 
ing the evolution of all flelds, especially by taking into 
account the back-reaction of the created particles on the 
inflaton oscillations and geometry. In any case, one sees 
that quantum fluctuations can be larger than (would be) 
thermal fiuctuations in equilibrium and thus they should 
be taken into account when the thermalization of the de- 
cay products is studied. 



Back- reaction 

Another interesting process that can potentially be 
affected by quantum fluctuations is the back-reaction 
of created particles on the evolution of the background 
fields. It is important to understand back-reaction effects 
to get a complete picture of the reheating process. Here, 
we consider the preheating model with interaction (37) 
reviewed above. When the created x particles are taken 
into account, the background scalar field equation in (32) 
should be modified as 



2 2 i 

g X (!) ■ 



0. 



(77) 



One usually invokes Hartree approximation and uses the 
expectation value < x^ > for x^ in (77). Since < x^ > 

docs not depend on spatial coordinates and since initially 
(i.e. just after the end of inflation) the (f) field is homoge- 
nous, only the zero mode continues to exist and one can 
ignore the spatial derivatives in (77). In that case, when 
the mean value of x^ grows to satisfy 

< >~ rn^, (78) 

the back-reaction effects become important. According 
to [15] this happens when 



rrr^ 



n-y ~ 



(79) 
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where riy^ is the average total number density of x parti- 
cles. 

However, as shown in the previous section, has fluc- 
tuations comparable to its vacuum expectation value, 
therefore as < > grows and (78) is satisfied, the "ac- 
tual value" of depends highly on the position, espe- 
cially when it is compared at scales larger than the corre- 
lation length £_c- To determine the subsequent evolution 
of the inflaton zero mode, wc imagine the space to be di- 
vided into regions of volume and analyze (77) in each 
region independently, trying to predict the full motion by 
gluing the results. Note that if only the zero mode sur- 
vives even when the fluctuations arc taken into account, 
then this should be a good approximation. 

From (77) and by ignoring the expansion of the uni- 
verse, the frequency of the oscillations in the i'th region 
is given by 



2 2,22 

W( .) = m +g Xi, 



(80) 



where x? denotes the value of x^ in that region. Since 
uji is (nc;arly) uniform in the i'th region, one can approx- 
imately write 



0j ~ $ sin {uJit) - 



(81) 



where denotes the restriction of (j) in that region. 

Let us now try to see if these local solutions, which are 
valid in different regions, can smoothly be glued. From 
(80), one sees that when g"^ < x^ >~ m^, the frequencies 
start to differ from one region to another as 



(5w 



1. 



(82) 



Therefore, in time f ~ 1/w ~ 1/m, nearly corresponding 

to a single average oscillation, the oscillations of the in- 
flaton field in different regions become completely out of 
phase. 

Assume now that (6 = J^^f/'i, where each (6, has its 
support in the z'th region. In that case, one can estimate 
the spatial derivatives in (77) as 



(83) 



By comparing with the mass term, the spatial derivatives 
can only be ignored if 



(e 



(84) 



It is not possible to satisfy (84) , when the inflaton passes 
through its minimum (/> = 0. When it reaches its maxi- 
mum = spatial derivatives can be neglected if 



m > 



1 
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(85) 



Using (45) as an estimate for ^'^^^ , (85) requires ^ 
1, which can only be satisfied if go ^ 10^, or so. As a 



result, one sees that when back- reaction effects become 
important the coherence of the inflaton oscillations is lost 
and the spatial derivatives must not be neglected. 

If preheating continues after back-reaction, our find- 
ings indicate a major change in the whole process. First 
of all, with the inclusion of spatial derivatives it becomes 
much more difficult, if not impossible, to determine the 
the subsequent evolution of the inflaton field. Moreover, 
one should also revise the particle creation process since 
the evolution of the background drastically changes. Pre- 
sumably, it will no longer be possible to talk about reso- 
nance bands in momentum space, but instead it would be 
more appropriate to analyze the process region by region. 



Formation of primordial black holes 

Finally, in this subsection we point out a new mecha- 
nism for the formation of primordial black holes, which 

can be efficient especially in models of preheating. As we 
will see, this process is different than the one considered, 
e.g. in [28], where a sufficiently large density contrast, 
whose amplitude is greater than a critical value as it en- 
ters the horizon, dynamically evolves in time to form a 
black hole. 

We show in the previous section that the energy den- 
sity p has fluctuations about its mean value. This is not 
surprising, since in quantum theory (4) implies that p be- 
comes a random variable as it is a function of a Gaussian 
random variable x- In general one can try to determine 
(at least an approximate) probability distribution for p, 
which must be close to a chi-squared distribution, but 
this is in general a difficult task since it contains two 
"non-commuting" random variables x and x- Instead of 
considering the full energy density, in the following we 
focus on the potential energy part 



1 



Pp 



2, ,2 



(86) 



which has a simple probability distribution determined 
in terms of the Gaussian distribution of x- It is clear 
that p > pp and it pp exceeds the critical value for a 
black hole to form then so does p. 

Consider now the potential energy Ep stored in a re- 
gion of size which is given by 



(87) 



If the Schwarzschild radius corresponding to this en- 
ergy becomes greater than then a black hole must 
form in that region (as discussed in the previous subsec- 
tions, £p^y' is in general smaller than the Hubble radius 
and thus it is legitimate to ignore the expansion of the 
universe). Recalling that the Schwarzschild radius for 
the energy E is given by Rs — E/M^, the condition 

Rs > ^,1^^" for the energy (87) becomes 



2 ^ 2 

X > Xo 



M2 



(88) 
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Since X is a Gaussian random variable with zero mean, 
the normahzed probabihty distribution can be written as 



nx) 



exp 



2<x2 > 



(89) 



Therefore, the probabihty of getting > Xo' which is 
also the probability for a black hole to form, is given by 



V = 



f-OO 

= 2 / P(X) 



dx- 



(90) 



Note that, we underestimate this probability, since only 
the contribution of the potential energy is considered. 
Nevertheless, (90) should give a good estimate since one 
would not expect a large hierarchy between the potential 
and total energies. 

On the other hand, because of the time dependence of 
the variance < > (and in some cases xo)j the proba- 
bility distribution (89) and the probability (90) changes 
in time. In that case, the black hole formation process 
can be visualized in time as follows. Assume at a fixed 
time t, some P{t) fraction of regions collapse to form 
black holes. Then one should wait for some time for the 
field X to evolve according to the new probability distri- 
bution and then regions containing large x^ fluctuations 
collapse again. Therefore, the process is actually cumula- 
tive. Moreover, the energy density in a collapsing region 
is larger than the average energy density in the universe, 
thus the actual fraction of energy that goes into primor- 
dial black holes must be larger than P. All these argu- 
ments support the use of (90) as a conservative estimate 
for black hole formation probability. 

If P denotes the observational upper limit of the frac- 
tion of energy that can go into the primordial black holes, 
then one should impose 



13 >V. 



(91) 



It is known that /3 < lO'^o (see, e.g. [28]). The 
magnitude of P depends very sensitively on the ratio 
< X^ > /Xo numerically one can find that to sat- 
isfy (91) with a (small) margin the following condition 
must be obeyed: 



X^ 



< X^ > 



> 22. 



Using the definition of Xo from (88), this implies 

/ 2 X K 
(< X >) < 



22M2(^r«^)2 



(92) 



(93) 



Therefore, in a model giving a large vacuum expectation 
value for x^ and a small correlation length S,^'^^^, for- 
mation of primordial black holes can be a problem. It 
is clear that this can be a dangerous issue especially for 
preheating. 

Let us check, for example, if the model (37) reviewed 
above passes the condition (93). Assume that the broad 



parametric resonance ends just when the back reaction ef- 
fects become important, i.e. q = 1 when < x^ >^ rn? . 
As shown in [15], for the realistic choice of the infla- 
ton mass m = 10~^Mp, this happens for the coupling 
5f ~ 3 X 10~^ and qq ~ 10^. Let us try to evaluate terms 
in (93) in terms of g. First of all, from (38), the mass 
parameter M can be seen to be oscillating with an am- 
plitude 5$. So, we take M = g$ in (93), i.e. we consider 
the black hole formation when the inflaton reaches its 
maximum value. As noted above, the vacuum expecta- 
tion value can be estimated as < x^ >— fn^ /g^ and the 
amplitude $ can be found from (7 ^ 1 as $ ~ "m/ g. Com- 
bining all these information, one sees that (93) implies 
g > 1.5 X 10~^, which is satisfied with a good margin since 
the actual value of the coupling is taken as 5 ~ 3 x 10"'^. 
Thus one concludes that formation of black holes is not 
an issue in this model. 



IV. CONCLUSIONS 

In this paper, we consider the well known particle 
creation effects in a time-dependent, homogenous and 
isotropic, classical background and point out a fea- 
ture that has not been elaborated in detail previously. 
Namely, we examine the fluctuations of important physi- 
cal quantities characterizing the particle creation process 
about their vacuum expectation values. We specifically 
consider the energy density, the pressure and the other 
components of the energy-momentum tensor, and find 
that all these quantities have in general large fluctuations 
comparable to average energy density. We also note that 
the deviation corresponding to the field square x^ also 
equals \/2 times the mean value of x^- It is possible to 
make sense of all these fluctuations by rising adiabatic 
or WKB regularization, therefore they have direct and 
unambiguous physical meaning. 

The spatial scale of these fluctuations is given by the 
correlation length of the quantum field x excited by the 
classical background. This is an important observation, 
which allows one to think of these fiuctuations being 
defined in regions that have the size of the correlation 
length, instead of imagining them as independent random 
variables defined at different points. If one still maintains 
the view that the field variables are defined point- wise, 
then one should always keep in mind that the variables 
inside the same region are correlated with each other and 
behave in the same way. 

In the second part of this paper, we focus on the re- 
heating process in single scalar field driven infiationary 
models and investigate possible implications of our find- 
ings. We consider two well known models, which are 
typical examples of the decay in broad parametric reso- 
nance regime and in perturbation theory, and determine 
the correlation length of quantum fluctuations. In both 
models, the physical correlation length becomes smaller 
than the Hubble radius at the end of the decay, where 
the ratio tends to be smaller in perturbation theory. 
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We investigate three possible effects of these fluctua- 
tions during reheating period. The first one is related to 
the final thernialization process of the decay products. 
One usually assumes that the decay products reach ther- 
mal equilibrium in a very short time, which is possible 
if the interaction rates are comparable to the expansion 
rate at the end of the decay. Although, this assump- 
tion is difficult to satisfy in realistic scenarios (see, e.g. 
[22]), it is important to estimate the maximum reheating 
temperature in a given model. We show that even when 
the interaction rates are presumed to be large enough 
and thermal equilibrium is expected to be set in a very 
short time, the existence of quantum fluctuations may 
delay this process. Namely, we observe that if on a given 
(subhorizon) scale quantum fluctuations are larger than 
would be thermal fluctuations in equilibrium then one 
should wait for energy to be redistributed in space to ob- 
tain the real equilibrium. Especially for larger temper- 
atures, the thermal correlation length can be many or- 
ders of magnitude smaller than the quantum correlation 
length, which requires more time for energy to spread. 
On the other hand, the whole process can actually be 
much more difficult to study since back-reaction effects 
must be considered to obtain the complete picture. 

As a second event which might be affected by quan- 
tum fluctuations, we consider the back-reaction of the 
created particles on the inflaton oscillations in the model 
of preheating. Back-reaction becomes important when 
the mean value of sufficiently grows to modify the fre- 
quency of the inflaton oscillations given by the inflaton 
mass. When this happens, because of the fluctuations in 
the oscillation frequency starts to appreciably vary on 
scales larger than the correlation length. We show that in 
a very short time the coherence of the oscillations is lost 
and spatial derivatives of the inflaton can no longer be 
neglected in determining its dynamical evolution. Prom 
that moment on, the whole back-reaction process changes 
and the problem becomes much more difficult to study 
since it involves non-homogenous fields in space. 

Finally, we notice that a large fluctuation in the energy 
density, which occurs in a region that has the size of the 
correlation length, causes a black hole to form if the cor- 
responding Schwarzschild radius becomes greater than 
the correlation length. To our knowledge, this process is 



considered for the first time in the context of reheating 
as a plausible mechanism for the formation of primordial 
black holes. In previous studies (sec e.g. [28]), the col- 
lapse of a sufficiently large, horizon size density contrast, 
which is known to produce a black hole as a result of its 
dynamical evolution under the influence of classical grav- 
ity, is considered as the main mechanism. We estimate 
the fraction of energy that can go into primordial black 
holes because of large energy fluctuations and show that 
observational constraints impose some new restrictions 
for the models. 

It is clear from our findings that, as far as the reheat- 
ing after inflation is concerned, the correlation length 
is one of the most important parameters characterizing 
the quantum fluctuations in the particle creation process. 
While sometimes a larger correlation length amplifies the 
impacts of quantum fluctuations as in the flnal thermal- 
ization process, in some instances a smaller ciorrelation 
length can increase the effects; formation of primordial 
black holes being an example as discussed below (93). 

It would be interesting to extend the present work in 
different directions. For example, it is of interest to study 
the back-reaction effects in more detail in the preheat- 
ing model. Particularly, new physics may arise if the 
broad parametric resonance continuous to exists when 
the back-reaction effects become important. It would 
also be interesting to extend the fluctuation analysis to 
include perturbations of the other flelds, i.e. the met- 
ric and the inflaton. The fluctuations in pressure and 
other components of the energy-momentum tensor are 
expected to play a role in this study. Finally, it is clear 
that the existence of these fluctuations should be con- 
sidered when symmetry breaking effects or formation of 
topological defects are studied. For example, in [29] for- 
mation of cosmic strings in a preheating model is studied 
by solving classical fleld equations numerically. Although 
in a linear theory such as considered in [29] , the classical 
field equations determine the evolution of the mean val- 
ues, it would not be surprising to see the destruction of 
cosmic strings longer than the correlation length because 
of the existence of large, order one fluctuations. It would 
be interesting to check this expectation by an explicit 
computation. 
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